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Abstract 



We use a numerical method to obtain the weak coupling perturbative coefficients 
of local operators with lattice regularization. Such a method allows us to extend 
the perturbative expansions obtained so far by analytical Feynman diagrams cal- 
culations. In SU (3) lattice gauge theory in four dimensions we compute the first 
eight coefficients of the expectation value of the Wilson loop on the elementary 
plaquette which is related to the gluon condensate. The computed eight coef- 
ficients grow with the order much faster than predicted by the presence of the 
infrared renormalon associated to the dimension of the gluon condensate. How- 
ever the renormalon behaviour for large order is quite well reproduced if one 
considers the expansion coefficients in a new coupling related to the lattice cou- 
pling by large perturbative corrections. This is expected since the lattice and 
continuum A scales differ by almost two orders of magnitude. 
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1 Introduction 



In asymptotically free theories, such as Yang-Mills gauge theories in four dimensions, the 
short distance distributions are reliably approximated by the perturbative expansion in 
at s (Q), the running coupling at the large momentum scale Q. However, even in this case the 
perturbative expansions are affected IQ-Q by ambiguities due to infrared (IR) renormalons, 
i.e. singularities on the integration contour in the Borel variable conjugate to a s (Q). The 
origin of these singularities is based on the renormalization group properties. Since in com- 
puting short distance distributions also small momenta k are involved in loop integrations, 
one has singularities due to the Landau pole in the running coupling a s {k) at k <C Q. 

The IR renormalons in the expectation value of composite operators have been studied 
in || in the case of the O(N) non-linear sigma model for iV — > oo. One finds that the related 
ambiguity is absorbed by non-perturbative contributions present in the various terms of the 
operator product expansion. A similar mechanism should also hold in the Yang-Mills case 
and prescriptions for resolving the ambiguities have been proposed || [| . 

Perturbative expansions of short distance distributions in QCD, such as the total e + e~ 
hadronic cross section and total decay widths, are known || up to three loops in the MS 
renormalization prescription. It is difficult to expect that one may reliably explore the 
properties of the IR renormalons on the basis of these not very long expansions. 

Recently we have proposed |7| a numerical method to obtain the weak coupling pertur- 
bative expansion of local observables in pure Yang-Mills gauge theory in four dimensions in 
the lattice regularization. The method is based on the idea of taking the weak coupling ex- 
pansion in the Langevin stochastic formulation of the theory M and solving numerically the 
truncated set of equations corresponding to a given perturbative order. By this method we 
hoped to be able to obtain perturbative expansions longer than the ones obtain by analytical 
methods @,g]. 

In this paper we present the results obtained in SU(3) lattice gauge theory in four di- 
mensions for the first eight coefficients of the expectation value of the "plaquette variable" , 
i.e. the Wilson loop on the elementary plaquette of lattice size a. The first four coefficients 
have been presented in Ref. |7j. On the basis of these new results we try to study the asymp- 
totic behaviour of the coefficients for large orders and we try to analyze the IR renormalon 
singularities. 

We consider the standard SU(3) gauge lattice action 



where a 5 is the usual coupling at the lattice scale. The sum extends to all plaquettes P 
in a hypercubic lattice in four dimensions, Up is the plaquette gauge field obtained from 
the link variable U^x) = exp{A fl (x) j V/3} . The observable we consider is the dimensionless 
expectation value of the plaquette 






Wi a \P) = 1 - | (Tr Up) . 



(2) 
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and we compute the coefficients c l £ l of the weak coupling expansion 

wt^HE^r", (3) 



71=1 



The physical interest of W±((3) is that, in the operator product expansion of this observ- 
able, the first non trivial term is the gluon condensate (a s F 2 ) which has dimension four. 
This quantity is renormalization group invariant and therefore its contribution to W&(0) is 
proportional to (aA) 4 , where A is the Yang-Mills scale given, at two loops, by 

(aA) 2 = C \JL\ e -^o f bo = 11/(47r) 2 f h = 102 /(47r) 4 . (4) 

There is no perturbative expansion for the quantity aA and therefore also for the physical 
condensate (a s F 2 ). The only perturbative contribution in (|3D is given by the first trivial term 
in the operator product expansion of W± (ft), i.e. the unit operator. The coefficients c 1 ® 1 are 
obtained from quartically divergent Feynman diagrams in which the lattice size a plays the 
role of the UV cutoff. The IR renormalons are present in this perturbative expansion since 
the momenta along the loops are integrated down to zero. 

The result of our calculation is that the first eight coefficients c 1 ^ grow much faster then 
expected from the conventional analysis of IR renormalons. For this comparison one uses 
the expansion 

Wl en {(5) = J2 c r n en f3~ n (5) 

71=1 

with c^ en determined by the presence of an IR renormalon singularity associated to the gluon 
condensate, i.e. of dimension four. We find then that c 1 ^ 1 grow with n much faster than c^ en 
for n < 8. This feature was already observed in the four loop calculation [[F] but now it is 
much more clear. Such a fast growth could be due to large subleading contributions in n. 
Indeed, as well known |It|, the lattice and continuum A scales differ by almost two order of 



magnitude so that the perturbative relation between the lattice and continuum couplings 

(3 C ont = Piatt -r + 0(P^ t ), (6) 

involves the large correction r = 12b Q ln(A/7rA^ ott ). This suggests to compare the lattice 
expansion Wl at (Piatt) with the renormalon expansion W£ en (j3 con t) in eq. @ where the two 
couplings are related by @. We find that for r ~ 2.4 the coefficients of in the expansion 
of Wl en (f3 con t) agree quite well with the coefficients cjf* for nk^A. Such a large value for r 
agrees with the values obtained flu]] in usual continuum renormalization schemes although 



it appears somehow larger than the one of MS renormalization. 

In Sect. 2 we briefly recall the method for obtaining numerically the perturbative coeffi- 
cients c^ GT for the SU(3) lattice gauge theory in four dimension. In Sect. 3 we present our 
results and discuss the statistical significance. In Sect. 4 we recall the asymptotic behaviour 
of the coefficients c^ GT deduced from the usual IR renormalon analysis by using the one 
and two loop running coupling. In Sect. 5 we discuss the behaviour of the eight computed 
coefficients in comparison with the one predicted by the IR renormalon analysis including 
the two loop correction of the running coupling. Sect. 6 contains some final comments. We 
include in the Appendix a conjecture, not plausible on the light of the present analysis, that 
the fast growth of the lattice coefficients c 1 ^ could be due, at least partially, to additional 
infrared singularities in the Feynman diagrams. 
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2 Weak coupling expansion on the computer 



Let us briefly recall the essential ingredients of our method for calculating the weak coupling 
expansion in the lattice regularization. First we describe the method for a scalar field; then 
we proceed to the case of SU (N) lattice gauge theory and its gauge fixing problem. (For a 
recent review about Langevin simulations see e.g. ||11||). 



The Langevin method consists in the simulation of a stochastic dynamical system having 
the field configuration space as its state space. Time evolution is dictated by the general 
equation 

where <fi is the field, S[(f>] the action and i] a Gaussian random noise satisfying the normal- 
ization 

(n(x, t)rj(x',t')\ = 5(x - x')S(t - t). (8) 

As a matter of fact, stochastic dynamics is devised in such a way that time averages along 
a trajectory converge to averages with respect to the Gibbs measure 

i £dt {Om)}) ri ^\\ D ^> Oft] e- 5M - (9) 
where the suffix rj denotes an average over the Gaussian noise. 

In order to implement Eq.^ on a computer, one can take t discrete with a time step dt = e 



<f>(x, n + 1) = <fi(x, n) — f(x, n), (10) 

where 

OS 

f(x,n)=e— r + y/ei](x,n) (11) 

0(f>{x, n) 

and now rj is normalized by: 

(i](x,n)r](x',n)) = d xx >8 nn >. (12) 

In this discrete form, Langevin equation is affected by a systematic error, which makes it 
necessary to extrapolate the results at e — > or to devise some higher order approximation 
to the continuum stochastic equation. In this paper we adhere to the simple recipe which 
consists in performing the simulation at several values of e with a linear fitting to e = 0. 

Let us now assume that the action is given by a free part (say So = ^J2 X 4>( X )(~^ + 
rrvfy(j){x)) plus an interaction term g Sx0 4 j the solution to Langevin's equation will then 
depend parametrically on the coupling constant g; we insert the formal expansion 

<f ) (x,t,g)~Y,9 k <P {k) (x,t) (13) 

fc>0 

into the Langevin evolution equation to obtain a hierarchy of stochastic evolution steps of 
the kind 

(f) {0) (x,n + l) = (f) {0) (x,n) -e(-A + ml)^ (x,n) + y/ei](x,n) 
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(p {1 \x,n + l) = (f) {1 \x,n) -e(-A + m 2 )(f) {1 \x,n) -s (4> (0 \x,n)) 

(2) (x,n+l) = <j) {2) (x,n) —e(-A + m 2 ,)^ (x, n)-e (0 (o) (x, nj) 2 (x, n) 



If we want to measure the expansion in g of any given observable 0[4>] we have just to 
insert the formal expansion (|H|) and take averages on the random process. For instance the 
expansion for the energy density will be given by 

(V0W) 2 ) + 2g • +g 2 ((V0«) 2 + 2V0<°> • V + ■ ■ ■ 



{£[<j>\) = ((V0(x)) 

The idea is very general and it can be applied to spin systems or to lattice gauge models. 

In £77(3) lattice gauge theory the action is given by Eq.(|I]) defined on a four dimensional 
periodic lattice. The link gauge variables U^(x) are SU(N) matrices labelled by the vector 
index \l = 1, . . .4 and by the site x of the lattice. Each configuration is then described by 
36 complex numbers at each lattice site. In this case the Langevin evolution step in discrete 
time e is given (see Ref Ji"3|) 



C^(z,n + 1) = e- F ^ x < n %(x,n), 



(14) 



where 



e(3_ 
AN 



x 



(15) 



and Hfj,(x) is a traceless antihermitian matrix extracted at random independently at each 
step from a Gaussian ensemble with normalization given by 



H ik (x,n)Hi m (x',n' 



H 



(16) 



In the weak coupling limit we can parameterize the link variables in terms of the usual gauge 
potentials A M , namely we set 



U ll {x) = e A "<>W, A^xY = -A^x), TiA, 
It follows that Langevin's equation takes the form: 



0. 



(17) 



e A„(x,n+l)/y/l3 



e -i^(x,n) e A^(a;,n)/v / /3 



where n is Langevin's discrete time. As in the scalar case the trajectories v4 M (x,n) will 
depend parametrically on (3: we expand them in a formal series in powers of /? -1 / 2 



J2P~ k/2 AP(x) 

k>0 



(19) 



and we rescale the time step e = t/(3 in such a way that Langevin's algorithm can be 
reformulated in terms of the fields Ajf- The formulae rapidly get rather cumbersome but 
they can easily be handled by a symbolic language or defining convenient structures as it is 
natural within the ZZ language of APE. 
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Figure 1: The eight independent time averages for each of the observables yielding an esti- 
mator of c[ GT , . . . , Cg GT at r = 0.01. 

We have reported the calculation of the average plaquette up to fourth loop (/3~ 4 ) in [0. 
It was shown that in order to obtain a stable process one has to impose some kind of gauge 
fixing; this fact seems quite natural in view of the infinite number of zero modes which one 
has to deal with in the expansion around a given classical vacuum. A gauge fixing procedure 



for Langevin quantization was proposed long ago by Zwanziger |J4| and it was implemented 
on the lattice in [[0|. The idea consists in introducing a new source in Langevin's equation 
which does not modify the asymptotic stationary distribution exp{— S} but it acts as an 
attractor towards the Landau gauge manifold defined by V ■ = 0. The implementation 
on the lattice consists in a gauge transformation which is executed after each Langevin step: 

U^x) -> e w[u ^(x) e ~ w[u ^ x+ ^\ (20) 

with 

w[U^x)] =\J2 (A-^ [U^x) - Ul(x - fi)])^ , A^U„(x) = U v (x) - U v (x - y), (21) 

A 4 



where the suffix "tr" denotes the traceless part of the matrix and A is a free parameter which 
we choose proportional to the time step r. 



3 Eight loop results for the plaquette 

The algorithmic setup presented in the last section lends itself to the possibility of performing 
high order weak coupling expansions; needless to say, since the calculation to n — th loop 
involves 2n — 1 auxiliary fields A^°\ A^\ . . . , A^ 2n ~ 2 \ it tends to be quite demanding in terms 
of computer memory and cpu time. The results reported here were obtained on an array 
processor with 64 units of the APE family. They required approximately 600 hours (roughly 
equivalent to 5000 hours on a single cpu of a Cray YMP) for the final run, which consisted 
of 24 independent Langevin histories of an average 4000 steps. Since the finite size effects 
on the single plaquette expectation value are rather small, we work on a 8 4 lattice with 
periodic boundary conditions which fits just fine on the memory of a single APE board, thus 
obtaining eight statistically independent Langevin histories at each run. There were three 
independent runs at r = .02, .015, .01 from which we extrapolate at the desired r = limit. 
We report in Fig. 1 the signal we get for r = 0.01. Since the higher loop observables depend 
on the lowest orders, we chose to start measuring the various observables in cascade, in order 
to let the system gradually thermalize. 

The statistical errors are computed by taking the standard deviation divided by square 
root of the number of statistically independent samples, which we estimate to be around 
64. This is a conservative estimate based on the presence of eight independent Langevin 
trajectories assuming an autocorrelation time of 300 steps. We plan to have a reliable mea- 
sure of autocorrelation in the near future, when longer Langevin histories will be available. 
Fig. 2 reports the process of extrapolation, together with the central value and the estimated 
statistical error for the first eight expansion coefficients. 



5 



n 


1 


2 


3 


4 


5 


6 


7 


8 


CLangevin 


1.998(1) 


1.218(1) 


2.940(5) 


9.28(2) 


34.0(2) 


134.9(9) 


563(5) 


2488(29) 


^analytic 


2 


1.218(7) 


2.9602 













Table 1: The expansion coefficients c„ obtained by the stochastic method and the known 
coefficients from Ref.(9). 

Figure 2: The extrapolation at r = of the first eight perturbative coefficients c^ GT on the 
8 4 lattice. 

We end this section by a remark about the finite size corrections. To get a rough estimate 
on them we performed a series of runs on a very small lattice (L 4 with L = 4) and on the 
largest feasible lattice on our APE (L = 12). We observe an increase in the values of c n (L) 
with L which is more sensible for higher loops. Essentially is appears rather reasonable that 
finite size corrections become important asn/L ~ 1. The overall picture seems to agree with 
the analysis of Ref. 0, but a deeper analysis would admittedly be desirable. 



4 IR renormalon analysis 

We recall here the usual conventional IR renormalon analysis. The Feynman diagrams of the 
plaquette expectation value W±{(3) are quadratically divergent. In general the perturbative 
contribution W2 a (f3) to the expectation of a composite operator which has UV divergences 
of degree 2cr (for the plaquette 2a = 4) has the form 

rQ 2 dk 2 / k 2 \ cr fi 

w ^=l t[<?) m - r 4 ™- (Q) ' (22) 

where Q is of the order of the inverse of lattice size, Q ~ 1/a, which plays the role of the 
UV cutoff. The function f(k) is dimensionless, renormalization group invariant and to first 
order is given by the one gluon exchange diagram, i.e. is proportional to a s (Q). In QED 
the contributions from the resummation of fermion loop diagrams |l]] gives f(k) = a(k), the 
running coupling in QED. The important simplification is that, from Ward identities, the 
vertex and wave function corrections cancel. In QCD the situation is different. Actually to 
reconstruct a s (k) one needs also vertex and wave function corrections which do not cancel. 
So in QCD one needs to consider more diagrams and there are no simple criteria to select 
the important ones. By taking into account that f(k) is renormalization group invariant one 
assumes that also in Yang-Mills gauge theory 

f(k) = a 5 (k) . (23) 

In this case the Landau pole in a s (k) is at a small value k -C Q and therefore the integral 
in (|22| ) diverges. The perturbative expansion of V^ 2 o-(/3) in ot s {Q) is obtained by expanding 
a s (k). The Landau pole disappears and each coefficient is finite. The divergence of (|22| ) at 
small k is reflected on the fact that the coefficients grow fast and the expansion diverges. A 
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nice framework for discussing this feature is given by the Borel transform. Introducing the 
(Borel) variable z 

■H 1 -!!)' (24) 

and using the integrand of ( p2|) can be written 

dk 2 fk 2 Y „, 4nb a 2 Jk) 6 , fl _ / 6 a 

where /3(a s ) is the beta function 

-l + 4x6„K W + ... (26) 



47r6 as(fc) 2 fro 

We have now to consider the range of integration for z corresponding to < k < Q. Beyond 
the Landau pole at z = z the mapping ( f24|) becomes ambiguous. However at two loops we 
have 

z = <Q) 1 1 + 4tt| a s (Q) | In ^ + • • • , 

thus we assume the range < z < oo. Neglecting the corrections in fl26|) we obtain the Borel 
representation 

with the perturbative coefficients given by 

WZTW = E < e "(2^ 61) /5" n , cr(2a, h) = C Ql±^ , 7 = (7 I , ^0 = ^ . 

n=i 2 o % o°o 

(28) 

where C is a numerical constant. A correction in the integrand in (|22[) proportional to a s (A;), 
as given for instance by taking into account the two loop contribution of the beta function in 
(p6f), simply modifies the normalization constant C . The IR renormalon singularity on the 
integration contour at z = z Q corresponds to the Landau pole of a s (k). For the plaquette we 
have 

2a = 4, ^ = ^~ = 4.785... (29) 
36 

To give a meaning to the representation (p7|) one has to give a prescription to specify the 
way the integration contour avoids the singularity z = zq. The various prescriptions differ 
by contributions proportional to the residue at the pole. The ambiguity in (]27 ) is then given 
by 

AW 8 (/3) ~ e"^ = e"^ ~ (aA) 2,J . (30) 

For the plaquette 2a = 4 this ambiguity is of the same order as the condensate ((a s F 2 ) ~ A 4 ), 
i.e. the first non trivial term in the operator product expansion of the plaquette. This sug- 
gests that the IR renormalon ambiguity should be absorbed by non perturbative contribu- 
tions coming from the various terms of the operator product expansion. 
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5 Renormalon and the numerical coefficients 



In Fig. 3 we plot the eight coefficients c 1 ^ of the perturbative expansion of the plaquette 
Wl at (f3) in d||) obtained from the numerical method reported in Sect. 3. We plot also the 
coefficients c r ^ n for the IR renormalon in (EH) with 2a = 4 and z is given in (P9"|). The 



normalization constant C is fixed in such a way that Cg en = c 8 a *. From the comparison we 
conclude that the growth of the coefficients c 1 ^ is much stronger than the one of c^ en (4, bi). 

A simple explanation for the discrepancy observed between c™ n and c^ t is that the 
coupling (3 = Pi a u in the lattice action ([I]) and the coupling suitable for the analysis of the 
IR renormalons in (|22"D are related by large perturbative corrections. One expects that the 
continuum coupling P con t is the suitable coupling for the renormalon analysis. As well known 
T0[| in the Yang-Mills case the A scales in the lattice and continuum renormalization schemes 



are quite different 

A cont = KAi a u, (31) 

with large values for K. In the MS scheme one has K = 28.8. If the coupling a s (Q) is 
defined by the static potential between quark-antiquark pair or by the three point vertex 
one finds the constants K = 46.1 or K = 69.4 respectively. This relation can be converted 
in terms of a relation between the lattice and continuum scale. From (U) and (31) one finds 



the relation (^J) between the two couplings where r = 1.85, r = 2.25 and r = 2.6 in the 
renormalization scheme in which the coupling is defined in MS, by the static potential or by 
the three point vertex respectively. 

In order to analyze the renormalon behaviour of the computed lattice coefficients c 1 ® 1 we 
compare the lattice expansion in the coupling Pf att with the renormalon expansion ( p8| ) 
in the coupling P~^ nt 



W^iPcont) = E C n (2a, 60 P~ n nt = E CT (2a, h) flS , (32) 

n=l n=l 

where we assume the following relation between the two couplings 



Pcont — Piatt — r — — , (33) 

Piatt 

so that C^ en (2a, b\) = c^ n (2a, b\) for r = r' = 0. The parameter r' takes contributions from 
three loops. We have r' = 6rb\/bo + 36<5&2/fro with 5&2 the difference between the three loop 
beta function coefficients in the lattice and continuum renormalization schemes. 

In Fig. 3 we plot the new coefficients C^ en (4, b\) obtained by fitting the two parameters r 
and r'. We find r = 2.41 and bbi = .12. For nk,A we have a very good agreement between 
c£*(4,&i) and C; en (4A). 

The value of r = 2.41 corresponds to the scale constant in Q3"Tp given by K — 56.1 which 
is between the ones obtained by defining a s {Q) by the static potential or by the three vector 
coupling and smaller than the one in the MS scheme. 

2 Notice that this value is remarkably close to Lepage and Makenzie's ([[l0|) taking into account the 
effective momentum q* ^2.8 proper to ay- 



S 



Figure 3: The values of c^*, (n = 1 . . . , 8) obtained by the stochastic method together with 
renormalon coefficients c^ en (4, bi) and c^ en (4, &i, r) with r = 2.41. 

Since the value of r is large, the relation between c 1 ® 1 and C^ en involves large cancelations. 
This implies that the results are sensible to r, as is rather evident from Fig 3. 

We have studied also the dependence on the two loop corrections by setting b\ = in the 
coefficients C^ en (4, b\ — 0). In order to obtain the fit with the lattice coefficients for n^4 we 
find r = 4.13 rather larger than the two loop value given before. 



6 Final considerations 



We used a numerical method to obtain on a computer long perturbative expansion in four 
dimensional Yang-Mills theory with a lattice regularization. The expansion parameter is the 
lattice coupling j3i aU = 2na s /3 entering in the Wilson action (]T|) at the lattice size a. The 
large order behaviour of the plaquette perturbative coefficients c^ ai agrees for n^4 with the 
presence of a IR renormalon at ( |2"9"D with 2cr = 4 corresponding to the dimension of the gluon 
condensate (a s F 2 ). However to obtain such a result one has to consider large corrections in 
the perturbative relation (|) between the lattice coupling (3i att and the continuum coupling 
f3 con t in which one discuss the renormalons. Such a large perturbative correction (r = 2.41) is 



expected since, as well known flO |, the A scales for the lattice and continuum regularizations 
differ almost by two orders of magnitude. It should be noticed that since c 1 ^ grow fast with 
n and r is large the coefficients in the continuum regularization are obtained in terms of the 
lattice ones by large cancelations. 

Since there are large perturbative corrections between lattice and continuum couplings 
it may be important to know the perturbative relation in (||) with many more coefficients. 
This could be attempted for instance by computing perturbatively by our method the static 
quark- ant iquark potential. 
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Appendix 

As discussed in the text the plausible explanation of the quite fast growth of the lat- 
tice coefficients c 1 ® 1 is the large perturbative corrections between the lattice and continuum 
couplings (||). However we would like to add an alternative explanation. The fast growth 
could be explained by an effective value of the dimension 2a much smaller than the actual 
value 2a = 4. Indeed for the effective value of 2a — 1 the coefficients c^ en (l,6i) in (|28|) fit 
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the numerical values of c 1 ^ 1 for n^3 — 4. To explain this smaller effective dimension suppose 
that, as suggested by perturbative QCD studies (see Jl6], [T7]]), the argument of the running 
coupling a s is modified by infrared corrections of Feynman diagrams in the Minkowski form 
so that the integrand f(k) is ( p3|) is given by 

f(k) = f 1 dxa 5 (xk) . (34) 
J o 

The perturbative contributions of f(k) in the expansion parameter a s (k) are finite and given 
by integrations over ln n x. Resumming the expansion in a s (k) one has a Landau pole in the 
variable xk instead of k. By assuming for the integrand in Q2%D the expression in ( |3"4"D one 
finds that the renormalon contribution in (1221) is 



(35) 

The first term, which gives the dominant contribution, is just the expression in (|22| ) with 
dimension 2a — > 1 with the Borel singularity at the universal value independent of the 
dimension a 

2a eff = l. z e Q ff = -4- = 1.196- •• (36) 
This Landau singularity would correspond to a perturbative ambiguity given by 

AW 2 {(3) ~ e - ^" ~ (aA) . (37) 
This large non perturbative corrections has been observed for quantities such as the average 
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The prove of the conjecture ([34]) would need a difficult analysis since one should study 
contributions in Feynman diagrams of terms with infrared singularities which are integrable. 
If the contributions in (|34]) are present one should care about the mechanism to resolve the 
perturbative ambiguity in (^). There are no operators to cure the ambiguity associated to 
dimension 2cr e // = 1. There is no contradiction with the analysis in the O(N) non-linear 
planar sigma model in two-dimensional H in which the position of the renormalon is actually 
related to the real dimension 2a. In two dimensions the structure of Feynman diagram 
singularities is completely different from the case of Yang-Mills theory in four dimensions 
especially if one uses the Minkowski form. The only connection between these models is the 
sign of the beta function, i.e. the ai"(Q) ln n k resummation. 
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